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ABSTRACT 

Pupil mapping is a promising and unconventional new method for high con- 
trast imaging being considered for terrestrial exoplanet searches. It employs two 
(or more) specially designed aspheric mirrors to create a high-contrast amplitude 
profile across the telescope pupil that does not appreciably attenuate amplitude. 
As such, it reaps significant benefits in light collecting efficiency and inner work- 
ing angle, both critical parameters for terrestrial planet detection. While much 
has been published on various aspects of pupil mapping systems, the problem of 
sensitivity to wavefront aberrations remains an open question. In this paper, we 
present an efficient method for computing the sensitivity of a pupil mapped sys- 
tem to Zernike aberrations. We then use this method to study the sensitivity of 
a particular pupil mapping system and compare it to the concentric-ring shaped 
pupil coronagraph. In particular, we show how contrast and inner working an- 
gle degrade with increasing Zernike order and rms amplitude, which has obvious 
ramifications for the stability requirements and overall design of a planet-finding 
observatory. 



Subject headings: Extrasolar planets, coronagraphy, Fresnel propagation, diffrac- 
tion analysis, point spread function, pupil mapping, apodization, PIAA 
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1. Introduction 

The impressive discoveries of large extrasolar planets over the past decade has inspired 
widespread interest in finding and directly imaging Earth-like planets in the habitable zones 
of nearby stars. In fact, NASA has plans to launch two space telescopes to accomplish this, 
the Terrestrial Planet Finder Coronagraph (TPF-C) and the Terrestrial Planet Finder In- 
terferometer (TPF-I) , while the European Space Agency is planning a similar multi-satellite 
mission called Darwin. These missions are currently in the concept study phase. In addition, 
numerous ground-based searches are proceeding using both coronagraphic and interferomet- 
ric approaches. 

Direct imaging of Earth-like extrasolar planets poses extremely challenging problems. 
To see why, consider viewing our solar system from a good distance (say from another 

nearby star system), our Sun would appear 10^° times brighter than Earth. Hence, we need 
an imaging system capable of detecting planets that are 10 orders of magnitude fainter than 
the star they orbit. Such a system is referred to as a high-contrast imager. Furthermore, 
given the distances involved, the angular separation for most targets is very small, requiring 
the largest launchable telescope possible. 

For TPF-C, for example, the current baseline design involves a traditional Lyot coro- 
nagraph consisting of a modern 8th-order occulting mask attached to the back end of a 
Ritchey-Chretien telescope having an 8m by 3.5m elliptical primary mirror (see, e.g., Kuch- 
ner et al. (2004)). Alternative innovative back-end designs still being considered include 
shaped pupils (see, e.g., Kasdin et al. (2003) and Vanderbei et al. (2004)), a visible nuller 
(see, e.g., Shao et al. (2004)) and pupil mapping (see, e.g., Guyon (2003) where this tech- 
nique is called phase-induced amplitude apodization or PIAA) . By pupil mapping we mean 
a system of two lenses, or mirrors, that take a flat input field at the entrance pupil and 
produce an output field that is amplitude modified but still fiat in phase (at least for on-axis 
sources) . 

Pupil mapping has received considerable attention recently because of its high through- 
put and small effective inner working angle (IWA). These benefits could potentially permit 
more observations over the mission lifetime, or conversely, a smaller and cheaper overall 
telescope. As a result, there have been numerous studies over the past year to examine 
the performance of pupil mapping systems. In particular, in Traub and Vanderbei (2003) 
and Vanderbei and Traub (2005), formulas for the optical surfaces were derived using ray 
optics. However, this analysis failed to account for the complete diffraction through a pupil 
mapping system. More recently, Vanderbei (2006) provided a detailed diffraction analysis. 
Unfortunately, this analysis showed that a pure pupil mapping system cannot achieve the 
needed 10~^° contrast; the diffraction effects from the pupil mapping systems themselves are 
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so detrimental that contrast attained at the first image plane is limited to 10~^. In Guyon 
et al. (2005) and Pluzhnik et al. (2006), a hybrid pupil mapping system was proposed that 
combines the pupil mapping mirrors with a slight apodization of oversized entrance and exit 
pupils. This combination docs indeed achieve the needed high-contrast point spread func- 
tion (PSF). In what follows, wc will refer to such apodized pupil mapping simply as pupil 
mapping, and use the term pure pupil mapping to refer to the original unapodized system. 

What remains to be answered is how pupil mapping behaves in the presence of opti- 
cal aberrations. It is essential that contrast be maintained during an observation while the 
system is being aberrated due to the small dynamic perturbations of the primary mirror. 
An understanding of this sensitivity is critical to the design of TPF-C or any other observa- 
tory. In Green et al. (2004), a detailed sensitivity analysis of both shaped pupils and Lyot 
coronagraphs is given. Kuchner et al. (2004) introduced an 8th order image plane mask for 
the Lyot coronagraph that reduces the sensitivity to low-order aberrations to the level of 
shaped pupils. Both of these design approaches achieve the needed sensitivity for a realizable 
mission. Up to now, however, no comparable study has been done for pupil mapping. 

In this paper, we present an efficient method for computing the effects of wavefront 
aberrations on pupil mapping. We begin with a brief review of the design of pupil mapping 
systems in Section 2. We then present in Section 3.1 a semi- analytical approach to computing 
the PSF of a pupil-mapping system in the presence of aberrations represented by Zcrnike 
polynomials. The assumption of Zernikes allows an analytical treatment of the azimuthal 
angle, reducing the computational problem from a double integral to a few single ones, thus 
eliminating the need for massive computing power. We also show how these methods can be 
used to compute PSFs for a purely apodizing coronagraph in Section 3.2 as well as off-axis 
PSFs in Section 3.3. We conclude Section 3 with a sensitivity analysis of pupil mapping 
systems in Section 3.4, showing the degradation of contrast and inner working angle with 
Zernike order and rms amplitude. We also demonstrate the achievable contrast limit of the 
current pupil mapping approach. In Section 4, we repeat the sensitivity analysis for an 
alternative coronagraph design, the shaped pupil coronagraph with concentric ring masks, 
and compare to pupil mapping. In particular, we show that the pupil mapping approach to 
high-contrast is more sensitive to aberrations than a purely apodized one, possibly making 
it more difficult to exploit the inner working angle advantage of pupil mapping. 
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2. Review of Pupil Mapping and Apodization 

2.1. Pure Pupil Mapping via Ray Optics 

We begin by summarizing the ray-optics description of pure pupil mapping. An on- 
axis ray entering the first pupil at radius r from the center is to be mapped to radius 
f = at the exit pupil. Optical elements at the two pupils ensure that the exit ray is 
parallel to the entering ray. The function R{r) is assumed to be positive and increasing or, 
sometimes, negative and decreasing. In either case, the function has an inverse that allows 
us to recapture r as a function of f: r = R{r). The purpose of pupil mapping is to create 
nontrivial amplitude profiles. An amplitude profile function A{r) specifies the ratio between 
the output amplitude at f to the input amplitude at r (in a pure pupil-mapping system the 
input amplitude is constant). Vanderbei and Traub (2005) showed that for any amplitude 
profile A{f) there is a pupil mapping function i?(f) that achieves this profile. Specifically, 
the pupil mapping is given by 



i?(f) =±y^ 2A'^{s)sds. (1) 

Furthermore, if we consider the case of a pair of lenses that are piano on their outward-facing 
surfaces, then the inward-facing surface profiles, h{r) and h{r), that are required to obtain 
the desired pupil mapping are given by the solutions to the following ordinary differential 
equations: 

dh , . r — R(r) , , , 

^(r) = ^ ' h{0)^z, (2) 

^JQl + in^-l){r-R{r)y 

and 

§(f)= . M0) = 0. (3) 

\J QI + (n" - ) - f)- 

Here, n is the refractive index and Qq is a constant determined by the distance z separating 
the centers (r = 0, f = 0) of the two lenses: Qq = —(n — l)z. 

Let S{r, r) denote the distance between a point on the first lens surface r units from the 
center and the corresponding point on the second lens surface f units from its center. Up to 
an additive constant, the optical path length of a ray that exits at radius f after entering at 
radius r = R{f) is given by 

Qo{f) = S{R{f),r) + n{h{f) - h{R{f))). (4) 

Vanderbei and Traub (2005) showed that, for an on-axis source, (5o(0 constant and equal 
to Qo- 
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2.2. High- Contrast Amplitude Profiles 

If we assume that a collimated beam with amphtude profile A{f) such as one obtains 
as the output of a pupil mapping system is passed into an ideal imaging system with focal 
length /, the electric field E{p) at the image plane is given by the Fourier transform of A{f): 

E{C,V) j j e-'"^'^A(VFT^)dydx. (5) 

Here, Eq is the input amplitude which, unless otherwise noted, we take to be unity. Since the 
optics are azimuthally symmetric, it is convenient to use polar coordinates. The amplitude 
profile A is a function of f = \Jx^ + y"^ and the image-plane electric field depends only on 
image-plane radius p = -^^2 _|_ ^2. 

E{p) = 1^ J I e-^'''^'''''^^-'^^A{f)fdedf (6) 

= ^ / ^0 (-2^^) Mr)rdr- (7) 
The point-spread function (PSF) is the square of the electric field: 

Fsi{p) = \E{p)\\ (8) 

For the purpose of terrestrial planet finding, it is important to construct an amplitude profile 
for which the PSF at small nonzero angles is ten orders of magnitude reduced from its value 
at zero. Vanderbei et al. (2003a) explains how these functions are computed as solutions to 
certain optimization problems. 

The high-contrast amplitude profile used in the rest of this paper is shown in Figure 1. 

We end this section by noting that the need to design customized amplitude profiles 
arise in many areas, usually in the context of apodization — i.e., profiles that only attenuate 
the beam. Slepian (1965) was perhaps the first to study this problem carefully. For some 
recent applications, the reader is refered to the following papers in the area of beam shaping: 
Carney and Gbur (1999), Goncharov et al. (2002) and Hoffnagle and Jefferson (2003). 



2.3. Apodized Pupil Mapping Systems 

Vanderbei (2006) showed that pure pupil mapping systems designed for contrast of 10~^° 
actually achieve much less than this due to harmful diffraction effects that are not captured 
by the simple ray tracing analysis outlined in the previous section. For most systems of 
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practical real- world interest (i.e., systems with apertures of a few inches and designed for 
visible light), contrast is limited to about 10^^. In Vanderbei (2006), certain hybrid designs 
were considered that improve on this level of performance but none of the hybrid designs 
presented there completely overcame this diffraction-induced contrast degradation. 

In this section, we describe an apodized pupil mapping system that is somewhat more 
complicated than the designs presented in Vanderbei (2006). This hybrid design, based on 
ideas proposed by Olivier Guyon and Eugene Pluzhnik (see Pluzhnik et al. (2006)), involves 
three additional components. They are 

1. a preapodizer Aq to soften the edge of the first lens/mirror so as to minimize diffraction 
effects caused by hard edges, 

2. a postapodizcr to smooth out low spatial frequency ripples produced by diffraction 
effects induced by the pupil mapping system itself, and 

3. a backend phase shifter to smooth out low spatial frequency ripples in phase. 

Note that the backend phase shifter can be built into the second lens/mirror. There are 
several choices for the preapodizer. For this paper, we use the preapodizer given by Eqs. (3) 
and (4) in Pluzhnik et al. (2006): 



where ^max denotes the maximum value of A{r) and /3 is a scalar parameter, which we take 
to be 0.1. It is easy to see that 



• Ao(r) approaches 1 as A{r) approaches Amax, and 

• Ao(r) approaches as A{r) approaches 0. 

Incorporating a post-apodizer introduces a degree of freedom that is lacking in a pure 
pupil mapping system. Namely, it is possible to design the pupil mapping system based on 
an arbitrary amplitude profile and then convert this profile to a high-contrast profile via an 
appropriate choice of backend apodizer. We have found that a simple Gaussian amplitude 
profile that approximately matches a high-contrast profile works very well. Specifically, we 
used 




• A{r)/A^^ < Ao{r) < 1 for all r. 



A 



pupmap 



(f)=3.35e-22(^/")'. 
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where a denotes the radius of the second lens/mirror. 

The backend apodization is computed by taking the actual output amphtude profile 
as computed by a careful diffraction analysis, smoothing it by convolution with a Gaussian 
distribution, and then apodizing according to the ratio of the desired high-contrast amplitude 
profile A{f) divided by the smoothed output profile. Of course, since a true apodization can 
never intensify a beam, this ratio must be further scaled down so that it is nowhere greater 
than unity. The Gaussian convolution kernel we used has mean zero and standard deviation 
a/ VlOO, 000. 

The backend phase modification is computed by a similar smoothing operation applied 
to the output phase profile. Of course, the smoothed output phase profile (measured in 
radians) must be converted to a surface profile (having units of length). This conversion 
requires us to assume a certain specific wavelength. As a consequence, the resulting design 
is correct only at one wavelength. The ability of the system to achieve high contrast degrades 
as one moves away from the design wavelength. 

Figure 2 shows plots characterizing the performance of the apodized pupil mapping 
system described in this section. The specifications for this system are as follows. The 
designed-for wavelength is 632. 8nm. The optical elements are assumed to be mirrors sepa- 
rated by 0.375m. The system is an on-axis system and we therefore make the non-physical 
assumption that the mirrors don't obstruct the beam. That is, the mirrors are invisible ex- 
cept when they are needed. The mirrors take as input a 0.025m on-axis beam and produce 
a 0.025m pupil-remapped exit beam. The second mirror is oversized by a factor of two; that 
is, its diameter is 0.050m. The postapodizer ensures that only the central half contributes to 
the exit beam. The first mirror is also oversized appropriately as shown in the upper-right 
subplot of Figure 2. 

After the second mirror, the exit beam is brought to a focus. The focal length is 
2.5m. The lower-right subplot in Figure 2 shows the ideal PSF (in black) together with 
the achieved PSF at three wavelengths: at 70% (green), 100% (blue), and 130% (red) of 
the design wavelength. At the design wavelength, the achieved PSF matches the ideal PSF 
almost exactly. Note that there is minor degradation at the other two wavelengths mostly 
at low spatial frequencies. 

It is important to note that the PSFs in Figure 2 correspond to a bright on-axis source 
(i.e., a star). Off-axis sources, such as faint planets, undergo two effects in a pupil mapping 
system that differ from the response of a conventional imaging system: an effective mag- 
nification and a distortion. These are explained in detail in Vanderbei and Traub (2005) 
and Traub and Vanderbei (2003). The magnification, in particular, is due to an effective 
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narrowing of the exit pupil as compared to the entrance pupil. It is this magnification that 
provides pupil mapped systems their smaller effective inner working angle. The techniques 
in Section 3 will allow us to compute the exact off-axis diffraction pattern of an apodized 
pupil mapped coronagraph and thus to see these effects. 

While the effective magnification of a pupil mapping system results in an inner working 
angle advantage of about a factor of two, it docs not produce high-quaity diffraction limited 
images of off-axis sources because of the distortion inherent in the system. Guyon (2003) 
proposed the following solution to this problem. He suggested using this system merely as a 
mechanism for concentrating (on-axis) starlight in an image plane. He then proposed that 
an occulter be placed in the image plane to remove the starlight. All other light, such as the 
distorted off-axis planet light, would be allowed to pass through the image plane. On the 
back side would be a second, identical pupil mapping system (with the apodizers removed), 
that would "umap" the off-axis beam and thus remove the distortions introduced by the first 
system (except for some beam walk — see Vanderbei and Traub (2005)). 

3. Diffraction Analysis of Apodized Pupil Mapping with an Aberrated 



A general 2D diffraction analysis of a pupil mapping system with a possibly aberrated 
input field requires significant computing power. In this section, we present a much more 
efficient method assuming the aberrations are Zernike polynomials, which allows us to treat 
the azimuthal variables analytically. In particular, we perform a full end-to-end diffraction 
analysis of an apodized pupil mapping system with aberrations, as well as with a tilted field 
such as one from an off-axis planet. Such an analysis is required to determine the sensitivity 
of these systems to phase errors at the input pupil, as well as the off-axis response. At the 
end of this section, we describe the results for specific cases of interest. 



If we assume that a flat, on-axis, electric field arrives at the entrance pupil, then the 
electric field at a particular point of the exit pupil can be well-approximated by superimposing 
the phase-shifted waves from each point across the entrance pupil (this is the well-known 
Huygens-Presnel principle — see, e.g.. Section 8.2 in Born and Wolf (1999)). For an apodized 
pupil mapping system, we can write this as 



Wavefront 



3.1. Mathematical Development 



E, 




XQ{x,y,x,y) 



1 



e^-'Qi^'y'^'y)/^Ai^{r)dydx, 



(9) 
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where 

Q{x, y, x,y) = ^{x- x)^ + {y - yf + {h{r) - h{f)Y + n{Z - h{r) + h{f)) (10) 

is the optical path length, Z is the distance between the piano lens surfaces (i.e., a constant 
slightly larger than z), Aya{r) denotes the input amplitude apodization at radius r, Aout(^) 
denotes the output amplitude apodization at radius f, and where, of course, we have used r 
and f as shorthands for the radii in the entrance and exit pupils, respectively. If the arriving 
field is not fiat but instead has a phase profile given by a{x, y), then the integral must include 
this phase shift: 

£^out(5^,^) = ^ut(r) / / -—A: -e2-«(™)/^+-(-'^)An(r)%cix. (11) 

J J ^wy^iy-i-^-iy) 

As before, it is convenient to work in polar coordinates: 

£;out(r,^~)=Ao,t(r) / / \ e^--'Q(->'^-^")/"+^"(-'^)An(r)rcigcir, (12) 

J J XQ[r,r,0 — 0) 

where 

g(f, r,e) = \lr'^- 2rfcos9 + P + {h{r) - h{f)Y + n{Z - h{r) + h{r)). (13) 

For numerical tractability, it is essential to make approximations so that the integral over 9 
can be carried out analytically. To this end, we need to make an appropriate approximation 
to the square root term: 

S = ^Jr'^ - 2rfcos9 + P + {h{r) - h{f)f. (14) 

We approximate the 1/Q{r,r,9 — 9) amplitude-reduction factor in Eq. (12) by the 
constant 1/Z (the paraxial approximation). The Q{r,r,9 — 6) appearing in the exponential 
must, on the other hand, be treated with care. Recall that Q{f, R{f), 0) is a constant. Since 
constant phase shifts are immaterial, we can subtract it from Q(f, r, 9) in Eq. (12) to get 

E,^,{f,9) ^ j j e2-(«(^>'^-^Ve(^.«(^)'0))A+-('-'^)A,(r)rd^dr. (15) 

Next, we write the difference in Q's as follows: 

Q{f,r,9-9)-Q{f,R{f),Q) = S{f,r,9-9)-S{f,R{f),Q)+n{h{R{f))-h{r)) 

S'^{r,r,9-9)-S'^{r,R{f),0) 



S(f,r,9-9) + S(f,R{f),0) 



+ n{h{R{f)) - h{r)ll6) 
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and then we expand out the numerator and cancel big terms that can be subtracted one 
from another to get 

S^{r,r,e -9) - S\r,R{r),0) = {r - R{f)){r + R{r)) - 2f (r cos{e - 6) - R{r)'^ 

+ {h{r) - h{R{r))) (h{r) + h{R{r)) - 2h{f)^ .(17) 

When r = R{f) and ^ = 0, the right-hand side clearly vanishes as it should. Furthermore, for 
r close to R{f) and 9 close to 9, the right-hand side gives an accurate formula for computing 
the deviation from zero. That is to say, the right-hand side is easy to program in such a 
manner as to avoid subtracting one large number from another, which is always the biggest 
danger in numerical computation. 

So far, everything is exact (except for the paraxial approximation). The only further 
approximation we make is to replace S{f,r,9 — 9) in the denominator of Eq. (16) with 
5'(f, R{f), 0) so that the denominator becomes just 2S'(f, i?(f), 0). Since S{r-, R{r), 0) appears 
many times in coming formulas, we abbreviate it as S{f) and hope this does not create any 
confusion. Putting this altogether, we get a new approximation, which we refer to as the 
Huygens approximation: 

i-'-'''^)l^^^-^^fi)A,^{r)d9rdr (18) 
j K{r,f)L{r,f,9)Ain{r)rdr, (19) 



Eou^{f,9) ^ / /e 



2m [ 

xe 



^out (^ 



xz 

where 

K{r, f) = e 



2m[ (^-fl(^))(^+HCT)+2^«(f)+(M^KM«(^)))(M^)+Mii(^))-2M^)) ) 
-2^ 



L{r,f,~9) = / %2.^^(-I^^^)/A+^aM)^^_ ^20) 

Jo 

Suppose now that the aberration is given by the (Z, m)-th Zernike polynomial 

a(r, 9) = eZl^ir/a) cos{m9), (21) 

where e is a small number. In this case, the kernel L{r,r,9) can be expressed in terms of 
Bessel functions. Recall that the n-th Bessel function Jn{x) is defined by the requirement 
that i"Jn{x) be the n-th Fourier coefficient in the Fourier series expansion 
is, 

Jn{x) = — ^ / e^^'^^^'^e*"^^^. (22) 
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The Fourier series then, of course, is simply 

oo 

gixcose ^ J2 i^Jk{x)e'^^. (23) 

k=—oo 

This series is usually referred to as the Jacobi-Anger expansion. 

Plugging the Zernike polynomial into the definition of the kernel L(r, f, 9) and substi- 
tuting a Jacobi-Anger expansion for each of the two cosine-exponentials, we get that 

L{r,fJ) = r%2.<-^^'"^)A+-M)rf^ (24) 
Jo 

= EE/ ^^'J^ ir^) ^''''-''MeZr{r/a))e^'-'d9 

3 k 

= J]e^^-^L,(r,f), (25) 

k 

where 

L,(r, f) = 27rz^(^-) (^^^ J, (6Zr(r/a)) (26) 
The second to the last equahty above follows from the trivial fact that 



/■27r 

Jo 



27r n = 

otherwise. 



Substituting Eq. (25) into Eq. (19), we can write 

= ^e^fe-^Eo„,,fc(f), (27) 



where 



Eo.t,k{f) = j K{r, f)L,{r, f)AUr)rdr. (28) 



Ignoring a complex unit exponential, the electric field in the image plane is just the 
Fourier transform of the electric field in the pupil plane. Working in polar coordinates, we 
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have 

Prom the definition (Eq. (22)) of the Bessel functions as certain Fourier coefficients, we get 

= 2ni''^e"'"''^Jkm(^-^rpy (30) 



Substituting Eq. (30) into Eq. (29), we get 

27r V ih^A, f ^ ( 27r 

k 



Y.^'^'^^E.^M. (31) 



where 



£^img,ik(p) ''""^'"^ ("0"''") (^2) 

Finally, since -Bimg,fe(p) = -E^img -fe(p), we see that 

oo 

-E^img(p,0) = -E'img,o(p) + 2^cos(A;m0)£;i„,g,fc(p). (33) 



A;=l 



Note that 

1 (X\k 



for < x <C 1. Hence, if we assume that e < 10~^, then the error incurred by dropping all 
terms in the sum on k for which A; > 2 will be approximately 10~^ at the most. Hence, any 
error in the PSF will be at the 10~^^ level. 

Thus, for a given Zernike polynomial aberration at the entrance pupil and e in Eq. (21), 
Eq. (33) can be used to compute the resulting response of the apodized pupil mapping 
system. Note that an arbitrary aberration can be expressed as a linear combination of 
Zernike polynomials and, for each Zernike, the azimuthal variable can be integrated explicitly, 
thereby reducing 2D integrals to sums of ID integrals. 
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3.2. Pure Apodization as a Special Case of Apodized Pupil Mapping 

If we set Apupmap(?^) = 1 for all f, then the pupil mapping system reduces to a trivial 
forward propagation of a flat beam. That is, the lenses/mirrors have perfectly fiat surfaces 
and no remapping occurs. Hence, in this case, the hybrid system becomes a pure apodization 
system with the desired high-contrast apodization being achieved by a pair of apodizers, one 
appearing before the beam is propagated and one after (except for the Presnel diffraction due 
to the propagation between the two mirrors.) In this fashion, the mathematics developed in 
3.1 can be applied to a pure pupil-apodization coronagraph. 

3.3. Off-Axis Diffraction Response 

In addition to modeling optical aberrations, the equations in Section 3.1 can be used to 
model off-axis sources such as planets because they are mathematically equivalent to large 
tilt aberrations. For example, a lau planet at lOpc is equivalent to a tilt aberration of 0.1 

arcseconds. For a D = 4 meter aperture and A = 632. 8nm, this would correspond to an 
angle of 3X/D and e = 9.63 in Eq. (21). Figure 3 is a contrast plot showing the result of 
computing Eq. (33) for such a planet, along with a simulation for a pure pupil-apodizing 
coronagraph. In addition, an on-axis star was added in both cases that is 10^ times brighter 
than the planet. We assumed a coronagraph entrance aperture of 25mm, and a focal length 
of 2.5m, so that for the apodizing coronagraph, or any conventional imaging system, the 
planet appears at roughly 194//m off-axis in the image plane. 

As mentioned before, pupil mapping both magnifies, which is good, and distorts, which 
is bad. Both effects are evident in Figure 3. The magnification moves the center of the 
distorted planet image out by roughly a factor of 2 (to about 400/im in Figure 3). Thus, 
the inner working angle of the system is effectively reduced, which is one of the benefits of 
pupil mapping. On the other hand, pupil mapping distorts the planet, which reduces the 
peak intensity a little with respect to the star's peak intensity. However, note that Figure 
3 shows contrast plots, i.e. intensity is normalized to peak star intensity for each case. On 
an absolute intensity scale, the whole pupil mapping curve will be higher because of the 
throughput advantage. Also note that the distortion caused by the two-mirror system can 
be countered by inverting the pupil mapping after starlight has been blocked, as described 
at the end of Section 2.3. 
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3.4. Analysis of Sensitivity to Zernike Aberrations 

Figures 4 to 7 show sensitivity results for the apodized pupil mapping system described 
in Section 2.3. Figure 4 is for aberrations that are 1/100-th wave rms whereas Figure 5 is for 
smaller 1/1000-th wave rms aberrations. Figure 6 shows radial cross-sections and Figure 7 
shows differences between the radial cross-sections of the aberrated and un-aberrated PSF's. 
(As noted, the (0, 0), or piston, Zernike is differenced with the ideal PSF.) Note that these 
plots show the PSF for an on-axis source, i.e. a star. By comparing these plots to Figure 3, 
one can directly compare the levels of the aberrations to the PSF of an off-axis planet. 

Two important effects are evident in Figure 7. First, the (0,0), or piston, Zernike 
represents the nominally achieved contrast and so its difference with the ideal shows the 
limitation of the current design approach. Any additional errors due to the Zernike aber- 
rations, as shown in the subsequent subplots (which are differenced with the nominal (0, 0) 
PSF), add to this static limit. Second, the pupil-mapping system shows a degradation in 
inner working angle with increasing Zernike and increasing aberration size. Any design that 
attempts to exploit the inherent magnification to achieve a smaller inner working angle must 
be able to achieve a higher level of aberration control. 

Another important issue to consider is the sensitivity to aberrations at wavelengths 
other than the design wavelength. Comparing Figure 6 with the bottom-right subplot of 
Figure 2, it is evident that the degradation due to aberrations at 1/1000 wave level are 
generally higher than the degradation due to detuning of the wavelength by a factor of 1.3 
or 0.7. Thus, the sensitivity to aberrations at these wavelengths should look very similar to 
Figure 6. Our simulations do indeed confirm this. Of course, at the longer wavelength, the 
aberrated PSF in Figure 6 would be expanded by the appropriate amount, and contracted 
for the shorter wavelength. 

4. Sensitivity Analysis for Concentric Ring Masks 

For comparison purposes, we now consider a high-contrast imaging system based on a 
simple concentric rings shaped pupil as described in Vanderbei et al. (2003b). A concentric- 
ring pupil-plane mask can be thought of as a circularly symmetric binary apodization. Hence, 
we assume that we are given such an apodization function A{r). 

As before, we assume that the input electric field could involve some aberration given 
by a{r, 9). The image plane electric field is given by the Fourier transform of the pupil- plane 
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field, which in polar coordinates is 



Expanding the two exponentials in a Jacobi- Anger expansion, we get that 



A{r)rdr. 



(35) 

The integral in brackets is easily solved, as before, eliminating one of the summations. With 
some rearranging, this yields 

i^img(p, = E ^('-™)'e^'=-^ Jkm (0rp) Jk {eZl\r/a)) A{r)rdT. (36) 

Finally, we use the symmetry property of the Bessel Functions {J-m{x) = {—l)"^Jm{x)) to 
yield the image plane field, 

oc 

-E'img(p, (p) = ^img,o(p) + 2 ^ cos{km(p)Ei^g^k{p), (37) 

k=l 

where 

^img,fc(p) = ^ ^t^'-'^^'Jkm (^^rp^ Jk {eZr{r/a)) A{r)rdr. (38) 
For a concentric ring mask having its edges at r^, the formula for £'img,A;(p) simplifies to 

E^r^M = E ^^'-"^^'Jkm (^rp) J, {eZr{r/a))rdr. (39) 

Figures 8 to 11 show sensitivity results for a high-contrast concentric-ring pupil mask. 
We also computed sensitivity results for the pure apodizing system; they appeared virtually 
identical to the results shown in these figures, which agree with the results in Green et al. 
(2004). Likewise, the off-axis results are indistinguishable from those of the pure apodizing 
system in Figure 3. 



Conclusions 



In this paper we presented an efficient method for calculating the distortions in the 
PSF of a pupil mapping system due to wavefront aberrations. Figures 4 to 11 show that our 
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particular pupil mapping system is somewhat more sensitive to low order aberrations than the 
concentric ring masks. That is, contrast and IWA degrade more rapidly with increasing rms 
level of the aberrations. However, this is partially mitigated by the magnification property 
of pupil mapping, making a direct comparison of the two systems more subtle, especially 
considering the distortion in pupil mapping. (A more direct comparison can be made for a 
pupil mapping that undoes the distortion, such as the unmapping system described at the 
end of Section 2.3, but that is outside the scope of this article.) In any case, it is evident that 
especially fine control of static and dynamic stability will be required in order to take full 
advantage of the smaller IWA intrinsic to the pupil mapping approach. A careful stability 
analysis of any design employing pupil mapping is thus necessary to determine its achievable 
operating range. 

Finally, we note that there is a spectrum of apodized pupil mapping systems. The two 
extremes, pure apodization and pure pupil mapping, both have serious drawbacks. On the 
one end, pure apodization loses almost an order of magnitude in throughput and suffers from 
an unpleasantly large IWA. At the other extreme, pure pupil mapping fails to achieve the 
required high contrast. There are several points along this spectrum that are superior to the 
end points. We have focused on just one such point, which was suggested by Guyon et al. 
(2005). We leave it to future work to determine if this is the best design point. With this 
paper, we have provided the tools to analyze the sensitivity of these kinds of designs. 
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Fig. 2. — Hybrid pupil mapping system as described in Section 2.3. Upper-left plot shows in 
red the target high-contrast amplitude profile and in blue the Gaussian profile that is a close 
approximation to it. The two black plots show the pre- and post-apodizers. Upper-right plot 
shows the lens profiles, blue for the first lens and red for the second. Lower-left plot shows 
in blue the phase map computed using Huygens propagation and in black the smoothed 
approximation. Lower-right plot shows in black the ideal PSF, in blue the actual PSF at the 
design wavelength, in red the actual PSF at 130% of the design wavelength, and in green 
the actual PSF at 70% of the design wavelength. (The horizontal scale for 130% and 70% 
has been dilated by factors of 1/1.3 and 1/0.7, respectively, and the x-axis labels correspond 
to the original A.) 
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Pupil Mapping 




Image-Plane Radius (mm) 

Fig. 3. — Distortion of off-axis sources in pupil mapping. Top: images of a 3X/D planet 
(e.g. a lau planet viewed from 10 parsecs away with a 4m telescope). The star appears as 
a saturated disk. Both images are contrast plots normalized to the peak intensity of the 
star image. The z scale is logio of contrast. Bottom: Horizontal slices. Note that these 
are contrast plots, i.e. everything is normalized to the peak intensity of the star image in 
pupil mapping as well as pure apodization. On an absolute intensity scale, the whole pupil 
mapping curve would be higher due to better throughput. 
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Fig. 4. — Pupil mapping. PSFs associated with various Zernikes normalized so that the 
rms error is, in each case, 1/100-th wave. The plots correspond to: Piston (0,0), Tilt (1, 1), 
Defocus (2,0), Astigmatism (2,2) Coma (3,1), Trefoil (3,3), Spherical Aberration (4,0), 
Astigmatism 2nd Order (4,2), Tetrafoil (4,4), etc. Important note: As explained in Section 
2.3, off-axis sources are "pushed out" by a magnification factor of about 2 in these high- 
contrast pupil-mapping systems. This needs to be remembered when comparing with later 
plots. 
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Fig. 5. — Pupil mapping. PSFs associated with various Zernikes normalized so that the rms 
error is, in each case, 1/1000-th wave. The plots correspond to: Piston (0,0), Tilt (1,1), 
Defocus (2,0), Astigmatism (2,2), Coma (3,1), Trefoil (3,3), Spherical Aberration (4,0), 
Astigmatism 2nd Order (4,2), Tetrafoil (4,4), etc. 
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Fig. 6. — Cross- sectional 1-D plots through the positive x-axis in the plots of the previous 
two figures. 
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Fig. 7. — Difference plots. For all the plots except (0,0), the difference is between the aber- 
rated and non-aberrated profile for pupil mapping. (The non-aberrated profile is essentially 
the (0,0), or piston term). For the (0,0) difference plot in the upper left, the difference is 
between the non-aberrated pupil mapping profile and the non-aberrated ideal apodization 
profile, showing that pupil mapping does not match the ideal apodization profile at the 
contrast level of 10^^^ and below. 
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Fig. 8. — Concentric Rings. PSFs associated with various Zernikes normalized so that the 
rms error is, in each case, 1/100-th wave. The plots correspond to: Piston (0,0), Tilt (1, 1), 
Defocus (2,0), Astigmatism (2,2) Coma (3,1), Trefoil (3,3), Spherical Aberration (4,0), 
Astigmatism 2nd Order (4,2), Tetrafoil (4,4), etc. 
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Fig. 9. — Concentric Rings. PSFs associated with various Zernikes normalized so that the 
rms error is, in each case, 1/1000-th wave. The plots correspond to: Piston (0, 0), Tilt (1, 1), 
Defocus (2,0), Astigmatism (2,2) Coma (3,1), Trefoil (3,3), Spherical Aberration (4,0), 
Astigmatism 2nd Order (4,2), Tetrafoil (4,4), etc. 
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Fig. 10. — Cross-sectional 1-D plots through the positive x-axis in the plots of the previous 
two figures. 
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Fig. 11. — Difference plots. Differences between the aberrated profile and the corresponding 
ideal. 



